Exact supersymmetry in the relativistic hydrogen atom in general dimensions 
— supercharge and the generalized Johnson-Lippmann operator 
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A Dirac particle in general dimensions moving in a 1/r potential is shown to have an exact J\f = 2 
supersymmetry, for which the two supercharge operators are obtained in terms of (a D-dimensional 
generalization of) the Johnson-Lippmann operator, an extension of the Runge-Lenz-Pauli vector that 
relativistically incorporates spin degrees of freedom. So the extra symmetry (S(2)) in the quantum 
Kepler problem, which determines the degeneracy of the levels, is so robust as to accommodate the 
relativistic case in arbitrary dimensions. 
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Introduction — There is an increasing amount of fas- 
cination with supersymmetry (SUSY) in various fields 
of physics. Q, In a broad context SUSY is a kind of 
(graded) Lie albebra that closes under a combination of 
commutation and anticommutation.Q One of the sim- 
plest realizations of such a symmetry may be found in one 
of the oldest problems in physics — motion of an elec- 
tron in a hydrogen atom, or, classically, Kepler's problem. 
It has long been known that there is a hidden symme- 
try (a dynamical symmetry) in the problem, which is 
related to an extra conservation (the Runge-Lenz vec- 
tor). So the symmetry of the problem is higher (SO (4)) 
than the trivial SO (3), which is the cause of the "acciden- 
tal" degeneracy (i.e., the energy level independent of I) in 
the spectrum of the hydrogen in non-relativistic quantum 
mechanics. If we move on to the Dirac electron in hydro- 
gen in the relativistic quantum mechanics, the accidental 
degeneracy is lifted (since the Runge-Lenz vector, which 
designates the direction of the perihelion, is no longer 
conserved). However, the degeneracy is lifted only incom- 
pletely, and two-fold degeneracies (for the Dirac-operator 
quantum number k = —(2S ■ L + 1) = ±(j + 1/2) with 
j being the total angular momentum) remain, which was 
puzzling. 

In 1985 Sukumar made an interesting suggestion that 
the strange degeneracy may be explained as a supersym- 
metry in the problem|j|. Subsequently however, Tanger- 
man et al.,[s| who have pointed out that there is indeed 
an exact AT — 2 supersymmetry for the nonrelativistic 
hydrogen atom, criticized that Sukumar's work has not 
actually constructed supercharges for the relativistic hy- 
drogen atom. On the other hand, analytic studies for 
the relativistic hydrogen in general spatial dimensions 
have been developed, and analytic solutions are now 
obtamedHQ- 

Given this background, the purpose of the present Let- 
ter is to show that an exact TV = 2 supersymmetry for 
the relativistic Dirac particle in a 1/r potential in fact 
exists in general D spatial dimensions. In doing so we 
have actually constructed the supercharges Q± (i.e., mu- 
tually aniicommuting operators that commute with the 



Hamiltonian) for the Dirac Hamiltonian in (D + 1) di- 
mensions. The symmetry enables us to obtain the low- 
est eigenenergy and its wave function for each sector of 
(the D-dimensional generalization of) \n\ in a simple and 
transparent manner, so the problem indeed turns out to 
be algebraically solvable. 

One interesting point is whether the supersymmetry in 
the hydrogen problem is related to the Runge-Lenz vector 
that describes the hidden conservation law on the non- 
relativistic level. Dahl et al.|j| have shown, for the rel- 
ativistic hydrogen in (3+1) dimensions, that the Runge- 
Lenz-Pauli vector (that necessarily involves spin degees 
of freedom for the Dirac particle, and is called Johnson- 
Lippmann operator 0) may be indeed used to construct 
supercharges describing the supersymmetry. So the ques- 
tion we address here amounts to whether this extends 
to general dimensions. We first construct the Johnson- 
Lippmann operator generalized to the relativistic hydro- 
gen in (D + 1) dimensions, and then show that the gener- 
alized operator may actually be used to construct the su- 
percharges. The supercharges are in fact shown to reduce 
to (the (D + 1) dimensional generalization of) Runge- 
Lenz-Pauli vector in the nonrelativistic limit. While the 
Runge-Lenz-Pauli vector in general dimensions has been 
discussed in Refs. ^flHl| for the nonrelativistic case, the 
relativistic one is obtained for the first time. So we shall 
conclude that the supersymmetry is unexpectedly robust 
enough to be extended to the most general case, i.e., rel- 
ativistic case in general dimensions. We can summarize 
the situation as 
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Dirac 's equation — Dirac's equation for the relativis- 
tic hydrogen atom in (D + 1) dimensions is written, in 
natural units (where c = h — 1), as 



{7°(p - eA Q ) + fPi ~ m}^(x) = 



(1) 



in standard notations, where 7's satisfy {7^, 7"} = g^ u — 



diag(l, -1, -1, ...), _p M = idfj 



„D 



), and 



2 



summations over repeated indices are implied. We as- 
sume a 1/r potential, so we have eAg = —Za/r, where 
a = e 2 /(Hc) is the fine-structure constant and Z the 
atomic number. When the lines of electric force are al- 
lowed to spread in the D spatial dimensions the Coulomb 
potential becomes l/r D ~ 2 , but we focus on the 1/r 
potential for general dimensions to retain the atomic 
stability 0|- We come back to this point later. If the 
time derivative is explicitly written, we have 



H = 7°m + 7°7V 



(2) 



Za 



(< = 1,2,. ..,D) 



where H is the Hamiltonian and summations are implied 
for repeated spatial superscripts. 

The operators that commute with the Hamiltonian are 
the total angular momentum (in D spatial dimensions), 



Jab = L a6 + - 7 V, (3) 
Lab = ix a db - ix b d a , 

and the Dirac operator (in D dimensions), 



K = 7 



(4) 



j°\j*-L 2 -S 2 + ^D-r 



that is related to the spin-orbit interaction. 

Generalized Johns on- Lippmann operator — While J 
and K are (usually the only) constants of motion for ar- 
bitrary central fields, it is known, for the ordinary D = 3, 
that there is an extra operator, a relativistic analogue 
of the Runge-Lenz-Pauli vector, that commutes with 
the Hamiltonian, as first constructed by Johnson and 
Lippmann[{| . We now generalize the Johnson-Lippmann 
operator to (D + 1) dimensions to define a Johnson- 
Lippmann-Katsura-Aoki operator, 



Zma 



A' 7 D+1 (i?- 7 V). (5) 



Here we have defined "f D+1 , a pseudo-scalar, which is 
a generalization of 7 5 in (3+1) dimensions, and which 
satisfies ( 7 D+1 ) t = 7 »+i,(y+i)2 = 1,{j d +\^} = 0. 
jD+i j g cons t ruc ted from { 7 °, 7 1 , . . . , r y D }, but its actual 
form depends on whether the spatial dimension is even 
or odd. 

From the (anti)commutation relations above, we can 
show, after a rather tedious manipulation, that we have 
indeed [H, A] = with a notable relation between H and 
A, 



K 



Za 



H 2 



(6) 



The expression reduces to the D = 3 counterpart ob- 
tained in Refs.|^.ll5|. 

For K and A, on the other hand, we can show that 

{A,K} = 0, 

i.e., we end up with mutually anticommuting operators, 
K, A, that commute with H . 

Construction of the J\f — 2 supercharge operators - 
We are now in position to construct the generators of the 
supersymmetry, which was left undone in 8] and stated 
111 Ref.@] as an operator interesting to find. We can do 
so by going back to an TV = 2 supersymmetric quantum 
mechanical model originally conceived by Witten[3T fl2T|. 
to which the present Hamiltonian is shown to be formally 
equivalent. Due to the relation ©, we can take H = 
A 2 as a Hamiltonian of the present problem. We can 
construct two operators Q\, Q2 that commute with 7i, 



AK 

Qi = A, Q 2 = i-p-p, 

\K\ 



(7) 



where k is the eigenvalue of K. We have then 

n = Qi = Ql (8) 
{Qi,Q 2 } = 0. 

If we further define 
Q± = \(Qi±iQ*) = Ul±^jA= ^A(l T V K ) (9) 



with V K = K/\k\, we have 



H 



= Q-=0, 
{Q+,Q-}- 



(10) 



This establishes an equivalence to Witten's model. 
So we can take the eigenvectors \n, ±) that satisfy 

H\n,±) = E&\n,±), (11) 
V K \n,±) = ±K±), 

since H= A 2 commutes with K, and we can talk about 
the simultaneous eigenvectors. 

V K does not commute with Q±, for which we have 
Q±\n, ±) = 0. On the other hand, [H, Q^] = im- 
plies that |n, ±) and Qz^\n, ±) have a degenerate eigen- 
value of H with different eigenvalues of V K . From the 
relation © a zero-eigenstate of H — A 2 is the ground 
state of the original H (or, more precisely, the lowest- 
energy state for each value of k). Since (n,a\TC\n,a) = 
(n,a\{Q+,Q-}\n,a) = \Q+\n,a)\ 2 + |Q_|n,a)| 2 > (for 
a = ±), a state is the ground state of H if the equality 
holds in the above inequality. The equality occurs when 
Q + |0, -) = or Q_|0, +) = 0. If we go back to the defi- 
nition of Q± (eq.@), the zero-eigenvalue state |0) should 
satisfy 

A\0) = 0, 



3 



since (l T p K )|n,T> ^0. 

Kernel of A - One way to establish the existence 
of such a state is an analytic method. Following Ref.[||, 
we write the wave function for odd spatial dimensions 
D = 2N + 1 as 



tp^x 1 ^ 2 , . . .,x D ) 



„_ N F(r)0 B (n) 
\iG(r)<P- K {n) 



(12) 



where F(G) is the "large (small)" components, <fi K is the 
angular part with angular coordinates f2, and Kip K = 
ki/j k holds. If we plug this into Aip K = 0, we have, after 
some manipulation, 



x— — h sign(n)x 




(13) 



= ^k 2 - (Zaf 



where we have defined F(r) = f(x), G(r) = g(x) in terms 
of a dimensionless x = (Zam/\n\)r ~ (Z/\K,\ao)r with clq 
being the Bohr radius. 

For this differential equation of second order, there 
are two independent solutions, f(x) oc x ~ s e~ B1Ba ^ x or 

xSe - S ign( K )x with 

K = ±[l+ 1 -(D-l)}, 

where l(— 0, 1, ...) is the orbital angular momentum, but 
the only normalizable one is the latter with k > 0. With 
a similar argument for g(x), we arrive at the kernel of A, 




FIG. 1: Energy level scheme for the relativistic hydrogen in D 
spatial dimensions. The lowest few levels are plotted against 
the Dirac quantum number k and tanhD, where the energy 
is plotted on a logarithmic scale to make the level splittings 
clearer. How the S(2) ladder is obtained with a shift in s and 
the supercharge operation Q+ is indicated by arrows with the 
sign of k and the principal quantum number n indicated. 



ip K oc x s N e 3 



(14) 



as the non-degenerate lowest state (the unpaired levels 
in Fig.l). 

For even D = 2N we can perform a similar proce- 
dure, again following Ref.0]. This time we can put 
ip K = r- N+1 l 2 [F(r)(j) K {n) + iG{r)cj>- K {n% for which 
the same differential equation results for F(r),G(r), 
so we have a non-degenerate ground state, ?p K oc 
x>-it+V*e-*{<f> K (n) + ^ ~ s)/Za]^ K {n)}. 

Eigenenergies and the group theory — Eigenenergies 
may readily be obtained algebraically: As used above, 
eq. © dictates that each of the zero-eigenstates of the su- 
persymmetric Hamiltonian TC — A 2 is the lowest-energy 
state (for each sector of K) of the original H . This im- 
mediately implies that E, the lowest-lying (within each 
sector of K) eigenvalue of H, is 



E/m 



Zct 

K 



Za 
s 



-1/2 



(15) 



in agreement with the analytical result in Ref. [6( . We can 
go up the ladder (where the leg corresponds to k = ±|«| 



while the rung spanned by a quantum number we call n') 
by making s = \J n 2 — (Za) 2 — ► s + 1 . 
So we end up with 



E/m 



(Zaf 



n - \k\ + + ^Jn 2 - (Za) 2 



-1/2 



where we assume Za < (D — 1)/2 for the atomic stability. 
Here the principal quantum number is given as 

n = I + 1 + n' 

with n' = 0, 1, .., where each level is doubly degenerated 
[corresponding to K = ±(Z + (D — l)/2) and related by 
Q+], except at the bottom of each ladder at n' — (i.e., 
n = I + 1) for which only k > should be taken (Fig.l). 
We can also see that the inter-dimensional degeneracy, 
which exists between the levels (7, D ) — > (I ± 1,D =p 2) 
noted for the nonrelativistic casejl3j, persists for the rel- 
ativistic case, and the supersymmetric ladders live on 
such a spectrum. 

Group-theoretically the present result implies the fol- 
lowing. The nonrelativistic hydrogen atom in D spatial 



4 



dimensions has a hidden symmetry (with the Runge-Lenz 
vector conserved), and the symmetry of the problem is 
higher (SO(L> + 1)) than the symmetry (SO(D)) of the 
space. If we go to the relativistic case the symmetry is 
degraded, but only partially degraded into SO(Z))(giS(2) 
due to the J\f = 2 supersymmetry. So the conjecture, 
stated in Ref.[]|, is established here for general dimen- 
sions. 

Relation with the Runge-Lenz-Pauli vector — Related 
to the above, we can note a certain relation between the 
Johnson-Lippmann-Katsura-Aoki operator (A) and the 
Runge-Lenz-Pauli vector as follows. Since A is Hermi- 
tian, we have 



A =I(A + At) = - 



(16) 



1 



2Zma 



X' 

r 



m r 



This expression reduces, in the nonrelativistic limit, to 

A -> -a l M\ 



where 



M l 



1 



ImZa. 



{p>U 



Ljip 1 ) 



is the nonrelativistic Runge-Lenz-Pauli vector in D spa- 
tial dimensions [ToL Hi) and a 1 = j d+1 j Y the spin oper- 
ator in D dimensions. Namely, the operator in this limit 
is the inner product of the Runge-Lenz-Pauli vector (M) 
and the spin (<x) in D dimensions flif. 

Summary and discussions — So the higher symmetry 
for the \jr potential in general dimensions is revealed 
to be surprisingly robust against the relativistic general- 
ization and is retained as a supersymmetry. The super- 
charge is indeed related to (a D-dimensional generaliza- 
tion of) the Johnson-Lippmann operator. 

Given the formula for the general D, we are tempted 
to ask what happens in the D — > oo limit. In the non- 
relativistic case the kinetic energy (~ (1/D)x potential 
energy) is dominated by the potential energy, so the sys- 
tem reduces to a set of harmonic oscillations around the 
classical potential minima as stressed in Ref . [l3| . In the 
relativistic case, however, we see that we cannot elimi- 
nate the small component (G) in the D — > oo limit, since 
the coupling between F and G does not vanish in this 
limit. This implies that we cannot trivially relate the 
supersymmetry with a set of oscillators even asymptoti- 
cally. 

Another basic question is whether the supersymmetry 
is just accidental to the 1/r potential. As mentioned 
above, electromagnetically there is a problem of how we 
can conceive the 1/r potential as the Coulomb potential 
in general dimensions. We should have l/r D ~ 2 potential 
from Gauss's law if the lines of electromagnetic force ex- 
tend over the D spatial dimensions, but there is a well- 
known Ehrenfest's 1920 result that the atom becomes 



unstable for this potential. Some authors argue that we 
should in fact stick to 1/r in general dimensions. ^} Ex- 
perimentally, a low-dimensional (D = 2) case may be 
interesting as an accessible one. 

In a broader context, a supercharge is a kind of gener- 
alization of the Dirac operator, so it could be that super- 
symmetry is shared by a wide class of Dirac-type equa- 
tions. Some authors [l^j have indeed discussed general- 
ization of the Runge-Lenz-Pauli vector, such as the case 
of a particle around a magnetic monopole with a vector 
potential around it. So the exploration of supersymmetry 
in wider gauge-field models will be an interesting future 
problem. We wish to thank Dr R. Arita for discussions. 
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